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ABSTRACT 

As  a  supplement  to  an  earlier  DSTO  Technical  Note,  Overview  of  STAP  algorithms,  this 
Technical  Note  examines  processing  of  phased  array  data  collected  from  non-side¬ 
looking  geometries.  Space-time  adaptive  processing  (STAP)  algorithms  developed  for 
the  radar  operating  in  the  side-looking  mode  cannot  in  general  be  directly  applied  to 
radar  operating  in  non-side-looking  modes,  including  the  forward-looking  mode.  A 
simple  Doppler  compensation  scheme  is  described,  so  that  the  Doppler  frequency  of 
mainlobe  clutter  is  shifted  to  zero  and  becomes  range  invariant,  like  the  clutter  data 
collected  from  the  side-looking  geometry.  After  Doppler  compensation,  all  STAP 
algorithms  developed  for  the  side-looking  geometry  can  be  applied.  Also  briefly 
discussed  are  some  techniques  for  suppression  of  range-ambiguous  clutter,  including 
so-called  3D  STAP  and  frequency  diversity  technology. 
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Doppler  Compensation  for  Airborne  Non-Side- 
Looking  Phased- Array  Radar 


Executive  Summary 

Australian  Defence  Force  (ADF)  has  been  equipped  with  the  Wedgetail  Airborne  Early 
Warning  and  Control  (AEW&C)  system  since  2009,  providing  an  entirely  new 
capability  for  surveillance  that  gathers  essential  and  necessary  information  from  air 
and  surface  surveillance  areas. 

One  of  key  sensing  elements  of  the  Wedgetail  AEW&C  system  is  the  airborne  L-band 
phased  array  radar  system.  When  its  side  arrays  are  operated  in  the  side-looking 
geometry,  the  earth  surface  has  a  Doppler  shift  which  is  range  invariant.  For  a 
homogeneous  clutter  environment,  the  collected  data  with  respect  to  range,  after  some 
normalisation  treatment,  can  be  considered  as  independent  and  identically  distributed 
(i.i.d.)  samples,  providing  a  means  for  estimating  the  covariance  matrix  of  unknown 
interferences.  The  covariance  matrix  or  the  like  is  essential  for  the  optimal  detection  of 
targets  of  interest.  Numerous  fully-adaptive  and  partially-adaptive  space-time 
adaptive  processing  (STAP)  algorithms  have  been  developed  for  the  processing  of  the 
phased-array  data  collected  from  this  side-looking  geometry. 

When  a  phased  array  radar  steers  away  from  the  side-looking  direction,  the  Doppler 
shift  of  mainlobe  clutter  becomes  range  dependent,  making  range  samples  non-i.i.d.  If 
a  STAP  algorithm  developed  for  the  side-looking  geometry  is  directly  applied  to  the 
data  collected  from  non-side-looking  geometries,  detection  loss  occurs  especially  in  the 
vicinity  of  the  mainlobe  clutter  Doppler.  In  addition,  there  are  also  more  false-alarms 
than  desirable  as  the  processor  cannot  sufficiently  suppress  the  mainlobe  clutter. 
Therefore,  the  STAP  algorithms  developed  for  the  side-looking  mode  cannot  in  general 
be  directly  applied  to  non-side-looking  geometries. 

Fortunately,  the  range-dependent  Doppler  shift  of  mainlobe  clutter  is  totally 
determined  by  the  parameters  of  radar  and  platform,  and  independent  of  the  clutter 
environment.  Following  previous  publications  in  the  literature,  this  Technical  Note 
describes  a  simple  Doppler  compensation  scheme.  Under  this  scheme,  data  collected 
from  a  non-side-looking  geometry  first  undergoes  a  Doppler  compensation.  After  the 
compensation,  the  Doppler  frequency  of  mainlobe  clutter  is  shifted  to  zero  and 
becomes  range  invariant,  like  the  data  collected  from  the  side-looking  geometry.  All 
the  STAP  algorithms  can  then  be  applied.  The  detection  improvement  by  this  Doppler 
compensation  scheme  has  been  demonstrated  using  the  data  generated  by  DSTO's 
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GPARM  (stands  for  Generic  Phase  Array  Radar  Modelling)  radar  system  simulation 
software. 

The  suppression  of  range-ambiguous  clutter  is  also  addressed.  Among  a  few 
techniques  discussed,  the  simplest  implementation  is  the  beamforming  in  elevation  in 
some  cases.  For  radar  looks  at  near  range,  depression  angles  of  the  mainlobe  clutter 
and  the  range-ambiguous  clutter  can  be  very  different.  A  beamforming  in  elevation,  if 
available,  can  null  the  illumination  of  the  range-ambiguous  clutter  to  achieve  clutter 
suppression.  3D  STAP  (azimuth,  elevation  and  Doppler  adaptive  processing)  can  be 
considered,  if  measurements  in  elevation  array  elements  are  available  and  the 
computational  cost  affordable.  The  3D  STAP  provides  extra  degrees  of  freedom  (DoF), 
so  the  clutter  coming  from  different  elevation  angles  can  be  suppressed.  The  drawback 
of  3D  STAP  is  the  significant  increase  in  the  requirements  of  sample  data  and 
computational  load. 

The  frequency  diversity  technology  where  different  channels  of  an  array  are 
transmitting  of  different  frequencies  but  on  receive  the  whole  array  receives  and 
processes  the  entire  transmitted  frequencies  has  the  potential  to  provide  suppression  of 
range  ambiguous  clutter.  A  frequency  diverse  waveform  generates  a  periodic  range 
dependent  beam  pattern  which  enables  range  ambiguous  clutter  to  be  managed  by 
careful  selection  of  the  range  ambiguity  interval  relative  to  the  underlying  periodicity 
of  the  range  dependent  beam  pattern  generated  by  the  frequency  diverse  waveform. 
Detailed  architectures,  characteristics  and  techniques  of  the  frequency  diverse  array 
radar  will  be  studied,  investigated  and  analysed  in  a  separate  report. 
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1.  Introduction 

Airborne  radar  systems  generally  seek  to  be  able  to  detect  targets  with  a  low  range  rate 
either  to  facilitate  the  detection  of  "beaming"  aircraft  or  for  the  detection  of  surface  moving 
targets  (i.e.  ground  moving  target  indication  (GMTI)).  For  airborne  phased  array  radar 
space-time  adaptive  processing  (STAP)  optimally  improves  the  detection  performance 
against  slow  moving  targets.  The  detection  of  slow  moving  targets  must  be  achieved 
against  a  background  of  radar  returns  from  the  Earth's  surface  (surface  clutter)  which  has 
an  apparent  Doppler  shift  due  to  the  interaction  of  the  finite  beam  width  from  the  antenna 
on  a  moving  aircraft.  The  scientific  literature  has  identified  a  number  of  adaptive  signal 
processing  techniques  to  produce  filters  that  will  minimise  the  impact  of  clutter  on  the 
radar  system,  and  many  of  these  techniques  require  training  data  from  independent  and 
identically  distributed  (i.i.d.)  samples  to  provide  the  statistical  information  to  optimise  the 
suppression  of  clutter  whilst  minimising  the  impact  on  potential  targets.  The  radar 
literature  (Klemm  2002;  2004;  Guerci  2003;  2010;  Ward  1994;  Melvin  2004;  Himed  et  al. 
2013)  primarily  considers  side-looking  radar  systems  because  for  a  radar  system  that  is 
looking  perpendicular  to  the  velocity  vector  the  dominant  clutter  return  produced  by  the 
antenna  mainlobe  intercepting  the  Earth's  surface  has  a  Doppler  spectrum  which  is  range 
invariant,  which  is  necessary  to  meet  the  requirement  for  "identically  distributed". 

The  operational  utility  of  an  airborne  radar  system  is  significantly  enhanced  if  the 
constraints  on  the  look  direction  of  the  radar  system  necessary  to  obtain  optimal 
performance  are  minimised;  in  this  case  the  ability  to  obtain  optimal  performance  when 
the  radar  is  scanned  well  away  from  the  side-looking  direction.  Some  of  the  literature 
(Wang  et  al.  1997;  Kreyenkamp  and  Klemm  2001;  Lapierre  et  al.  2009)  has  examined  non¬ 
side-looking  geometries.  In  this  report  we  further  develop  these  techniques  to  extend  the 
performance  of  algorithms  developed  for  the  side-looking  case  to  a  wider  angular 
coverage. 

All  fully-adaptive  and  partially-adaptive  STAP  algorithms  inclusively  need  to  estimate  the 
covariance  matrix  of  unknown  interferences  or  the  like  to  implement  the  adaptive 
processing.  Ideally  the  estimation  requires  range  samples  to  be  i.i.d.  Unfortunately,  the 
data  collected  from  a  non-side-looking  geometry  are  Doppler  range-dependent  and  hence 
not  i.i.d.  even  if  the  clutter  environment  itself  is  homogenous.  Use  of  non-i.i.d.  samples  for 
estimating  the  covariance  matrix  or  the  like  violates  the  estimation  requirement. 

The  fundamental  difference  between  the  clutter  data  collected  from  non-side-looking  and 
side-looking  geometries  lies  in  that  the  former  in  general  exhibits  range  dependent 
Doppler,  violating  the  i.i.d.  assumption  that  is  implicitly  applied  to  techniques  and 
algorithms  developed  for  the  side-looking  geometry.  Fortunately,  the  Doppler  response  of 
the  Earth  surface  clutter  is  determined  by  the  radar  and  platform  parameters  making  it 
possible  to  compensate  the  Doppler  from  range  bin  to  range  bin  in  such  a  way  that  the 
Doppler  appears  to  be  statistically  identical  for  all  range  samples,  like  the  Doppler  range 
invariant  data  collected  from  the  side-looking  geometry.  Once  the  Doppler  compensation 
is  done,  the  data  can  be  heated  as  if  they  were  collected  from  the  side-looking  geometry. 
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and  hence  all  the  techniques  and  algorithms  developed  for  that  geometry  can  be  applied 
directly. 

It  is  worth  noting  that  the  role  of  GMTI  may  not  be  compromised  for  the  non-side-looking 
mode1.  Although  the  mainlobe  clutter  has  a  nonzero  Doppler  shift  for  radar  looking  into  a 
non-side-looking  direction,  it  is  only  the  relative  velocity  of  target  to  the  ground  matters. 
The  relative  motion  between  the  radar's  main  beam  and  the  ground  only  results  in  a 
Doppler  shift  but  does  not  change  the  Doppler  difference  between  the  target  and  the 
clutter,  assuming  the  target's  radial  velocity  relative  to  the  ground  to  remain  the  same. 

This  Technical  Note  is  a  supplement  of  an  earlier  Technical  Note  (Dong  2011),  to  include 
the  processing  of  non-side-looking  data  for  AEW  radar.  Section  2  describes  the  range- 
dependent  Doppler  of  non-side-looking  data  and  its  compensation.  Section  3  presents  a 
numerical  example  showing  the  effectiveness  of  the  Doppler  compensation  scheme 
proposed  in  this  note.  In  addition,  range-ambiguous  cases  are  discussed  in  Section  4.  Also 
discussed  in  Section  5  are  techniques  of  3D  STAP  and  the  frequency  diversity  technology 
aiming  at  suppression/  mitigation  of  range-ambiguous  clutter  to  improve  the  performance 
of  GMTI.  Section  6  concludes  the  note. 


2.  Range-Dependent  Doppler  and  Its  Compensation 

Figure  1  depicts  geometries  of  forward-looking  and  side-looking  airborne  arrays.  The 
ground  is  approximated  to  a  horizontal  plane  and  the  platform  flies  horizontally  in  the  x 
direction  with  a  velocity  of  va .  The  forward-looking  array  is  the  array  whose  boresight  is 
in  the  x  direction  and  the  side-looking  array  is  the  one  with  its  boresight  in  the  y 
direction.  Owing  to  the  relative  movement  between  the  array  and  a  stationary  scatterer 
located  at  P  on  the  ground,  the  Doppler  frequency  of  the  scatterer  measured  by  the  radar 
is  (Klemm  2002,  Chapter  1;  Ward  1994), 


fa  =-£-™e™v-*a) 


(1) 


where  A  is  the  wavelength  of  carrier  frequency,  9  is  the  depression  angle  (under  the 
assumption,  it  is  also  the  grazing  angle),  (f>  is  the  azimuth  angle  and  <t>a  is  the  aircraft's 
drifting  (crabbing)  angle,  i.e.,  the  actual  flying  direction  is  away  from  the  x-direction  by  (j>a 
normally  due  to  the  atmospheric  wind  effects. 


1  This  is  especially  true  for  radar  operated  in  the  range-unambiguous  low  PRF  mode. 
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Figure  1:  Geometry  of  forward-looking  and  side-looking  airborne  arrays. 


The  depression  angle  is  a  function  of  the  ratio  of  platform  height  h  to  range  r ,  resulting 
in. 


(2) 


where  fa  niax  =  2vfl/2  is  the  maximum  Doppler  frequency  the  ground  clutter  can  have.  It 
can  be  seen  that  the  Doppler  frequency  is  dependent  on  the  ratio  of  hlr  whenever 
cos(fZ J  -  </y )  *  0  (i.e.  non  side-looking  situations). 

The  normalised  surface  clutter  Doppler  (isodops)  of  the  whole  surface  is  shown  in  Figure  2 
where  the  platform  flies  horizontally  in  the  x  direction  without  any  crabbing  angle  (</>a  =  0 ), 
coordinates  of  (x,  y)  are  also  normalised  by  the  platform  height  h  .  The  maximum  positive 
Doppler  and  negative  Doppler  are  reached  at  (oo,0)  and  (-oo,0) ,  respectively,  whereas  zero 
Doppler  is  kept  along  the  whole  ±y  directions  (the  side-looking  geometry). 

Shown  in  Figure  3  is  the  range-dependency  of  normalised  Doppler  for  mainlobe  clutter 
when  the  radar  steers  towards  three  non-side-looking  directions,  namely,  0°  from  nose 
(forward-looking),  30°  from  nose  and  60°  from  nose,  respectively.  It  can  be  seen  that  the 
range-dependent  Doppler  varies  significantly  at  near  ranges  (low  r/h  ratios).  At  far  ranges 
(e.g.  r/h>  10),  the  Doppler  of  mainlobe  clutter  becomes  approximately  range  invariant 
although  the  magnitude  of  the  Doppler  is  still  dependent  on  the  look  direction. 
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Figure  2:  Normalised  surface  clutter  Doppler  (the  isodops)  when  the  platform  flies  horizontally  in 
the  x-direction. 


r/h 


Figure  3:  Range  dependency  of  the  normalised  Doppler  frequency  ofmainlobe  clutter  (0°  from  nose 
is  the  case  of  radar  looking  in  the  flight  direction). 
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Shown  in  Figure  4  are  two  half  surface  rings,  corresponding  to  two  different  ranges,  for 
radar  operated  in  the  forward-looking  geometry.  Clutter  signal  received  by  antenna 
element  « ,  pulse  m  at  range  r  is  an  integral  of  the  echoes  from  the  constant  range  ring 
(returns  from  ambiguous  range  are  not  included  here  and  will  be  discussed  later),  given 
by. 


Zr(n,m) 


fa=n!2 

n  j.  \  l2.Kn  /„  i2nm  tnTr 

Z  $(0,<j)k)eJ  JseJ 

fa  =—7T  /  2 


for  «  =  0, 1, •  ■  • , TV  - 1  and  m  =  0,l,---,M -1 


(3) 


where  ^(6,<j>k)  is  the  backscattered  field  from  an  clutter  patch  area  of  A Ak  =  rcos8ArA</>k  at 
(0,<pk)  (range  r  and  height  h  determine  the  depression  angle  Q),  A r  is  the  radar's  range 
resolution;  Tr  =  1/  fPRF  is  the  pulse  repetition  interval  (PRI)  and  fPRF  is  the  pulse  repetition 
frequency  (PRF);  fs  and  fa  are  the  spatial  and  temporal  frequencies,  respectively.  The 
reason  for  the  integration  of  a  half  ring  only  is  that  the  backlobes  of  the  antenna  are 
assumed  to  be  small  enough  to  be  ignored.  Quantities  in  (3)  are  further  discussed  in  the 
following. 


Figure  4:  Clutter  signal  collected  by  an  antenna  element  for  a  pulse  is  an  integral  of  scattering  fields 
from  all  scatterers  in  a  half  ring  (tivo  rings  corresponding  to  two  different  ranges  are 
shown). 


The  backscattered  field  may  be  written  as  ^(0,</>k)=\^6,</>k)\ej,pk  where  q>k  is  the  initial 
phase  of  the  clutter  patch  relative  to  the  reference.  Since  phase  centres  of  scatterers  are 
random,  the  initial  phase  can  be  assumed  to  be  random,  if  clutter  patches  are  mutually 
uncorrelated.  The  amplitude  of  backscattered  field  may  be  determined  from  the  radar 
equation  should  the  backscattering  coefficient  be  known,  as  (Ward  1994), 
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P,TpG(0, </>k  )g(6,  <K  )23<t„  (A  <j>k  )AAk 
(4  n?Lsr4 

where  Pt  is  the  transmitted  power,  Tp  is  the  pulse  width;  G  and  g  are  the  full  array 
transmit  gain  and  the  element  receive  gain,  respectively;  <r0  is  the  backscattering 
coefficient  and  Ls  denotes  the  total  loss  for  the  radar  system.  The  single-element  single¬ 
pulse  clutter-to-noise  ratio  (CNR)  is, 

CNR  =  e{\%\2}/ct2  (5) 


where  a2  is  the  variance  of  thermal  noise. 

For  an  array  operated  in  the  forward-looking  mode  and  flying  in  the  x-direction  with  no 
crab  angle  ( <j)a  =  0 ),  the  spatial  and  temporal  frequencies  are,  respectively  (Klemm  2002; 
Ward  1994), 

fa  =^f-COS0COS<l>k  (6) 

fs  =^e-es)^k  (7) 

A 

where  9S  is  the  array  tilted  angle.  Often  radar  is  pointed  to  the  area  of  interest  for 
detection,  which  means  its  mainlobe  is  pointed  to  the  area  of  interest,  so  that  cos(<9 - 9S)  ~  1 . 

As  indicated  by  (6),  the  Doppler  frequency  of  mainlobe  clutter  ( cos  <j>k  =1),  when  the  radar 
system  is  operated  in  the  forward-looking  mode,  is. 


2v„  „  2v  1 


/a0=-^co  S0  =  -S1- 


2 


2 


h 

yrj 


(8) 


Equation  (8)  clearly  shows  the  range-dependency  of  the  mainlobe  clutter  for  radar 
operated  in  the  forward-looking  mode.  While  the  mainlobe  clutter  shall  be  dominant  in 
the  received  clutter  signal,  contributions  of  clutter  returns  from  other  azimuths  should  not 
be  ignored,  and  their  relative  strengths  are  largely  determined  by  transmit  and  receive 
gains  as  well  as  the  backscattering  coefficient  as  indicated  by  (4). 

For  the  purpose  of  demonstration.  Table  1  lists  parameters  of  a  simple  airborne  linear 
phased  array  used  in  the  simulation.  According  to  parameters,  the  transmit  beam  pattern 
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(the  whole  array  beam  pattern)  and  the  element  receive  beam  pattern  are  shown  in  Figure 
5  and  Figure  6,  respectively. 


Table  1:  Parameters  used  in  the  simulation 


Parameter 

Value 

Platform  height,  h 

5000  m 

Platform  velocity,  va 

150  m/s 

Pulse  repetition  frequency,  fPRF 

3200  Hz 

Carrier  waveform  length,  A0 

0.25  m 

Bandwidth  of  LFM 

5  MHz 

Sampling  rate 

10  MHz 

Number  of  array  elements,  N 

25 

Number  of  pulses  in  a  CPI,  M 

31 

Antenna  element  beam  pattern  in  azimuth 

cos  <p 

Figure  5:  One-ivay  full  array  azimuth  beam  pattern  (transmit)  used  in  the  simulation. 


Figure  6:  One-way  element  azimuth  beam  pattern  (receive)  used  in  the  simulation. 
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Using  the  parameters  given  in  Table  1,  clutter  received  by  receive  elements  are  simulated 
using  (3).  In  the  simulation,  the  dependency  of  clutter  intensity  upon  range,  as  shown  in 
(4),  is  normalised.  The  backscattered  field  is  determined  assuming  the  backscattering 
coefficient  cr0  to  be  a  constant  and  the  associated  initial  phase  to  be  random  but  uniformly 
distributed  in  [0,  2k]  .  Therefore,  the  clutter  intensity  received  from  ring  r,  and  ring  r2 
shall  be  statistically  identical  with  Doppler  frequencies  being  different  and  dependent  on 
range  as  shown  in  (6).  The  single-element  single-pulse  CNR  is  set  to  30  dB.  The  simulated 
data  for  each  range  bin  are  stored  in  a  two-dimensional  NxM  matrix,  where  N  and  M 
denote  the  number  of  array  elements  and  the  number  of  pulses  in  a  coherent  processing 
interval  (CPI),  respectively. 

Figure  7  shows  the  pulse-based  M  beamformed  clutter  spatial  spectra  (i.e.,  each  column  of 
the  NxM  matrix  is  beamformed)  using  the  simulated  clutter  data  for  a  radar  operating  in 
the  forward-looking  mode  with  r/h  =  2.0  (a  near  range  case).  The  channel-based  N 
Doppler  spectra  (i.e.,  each  row  of  NxM  matrix  is  Fourier  transformed)  are  shown  in 
Figure  8.  For  comparison,  the  spatial  spectra  and  Doppler  spectra  for  r/h  =  20  (a  far  range 
case)  are  shown  in  Figure  9  and  Figure  10,  respectively.  It  can  be  seen  that  the  spatial 
spectra  are  statistically  identical  for  the  near  range  and  far  range  cases.  The  only  difference 
between  the  two  cases  is  their  Doppler  spectra  that  peak  at  1039.2  Hz  for  r/h  =  2.0  and 
1198.5  Hz  for  r/h  =  20,  respectively,  according  to  (8).  It  is  worth  noting  that  in  the 
simulation,  the  backscattering  coefficient  cr0  is  assumed  to  be  a  constant  and  the  associated 
phase  to  be  random  with  a  uniform  distribution  in  [0,  2k]  . 


Figure  7:  Pulse-based  M  spatial  clutter  spectra  beamformed  from  the  simulated  clutter  data  for  a 
radar  operating  in  the  forward-looking  mode  with  r/h  =  2.0  . 
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Figure  8:  Channel-based  N  Doppler  spectra  of  the  simulated  clutter  data  for  a  radar  operating  in  the 
forward-looking  mode  with  r/h  =  2.0  (The  Doppler  frequency  peaked  at  1039.2  Hz). 
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Figure  9:  Pulse-based  M  spatial  clutter  spectra  beamformed  from  the  simulated  clutter  data  for  a 
radar  operating  in  the  forward-looking  mode  with  r/h  =  20 . 
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Figure  10:  Channel-based  N  Doppler  spectra  of  the  simulated  clutter  data  for  a  radar  operating  the 
in  forward-looking  mode  with  r!h  =  20  (the  Doppler  frequency  peaked  at  1198.5  Hz). 


If  the  Doppler  spectra  shown  in  Figure  8  is  circularly  shifted  by  a  proper  amount,  so  that 
their  peak  have  the  same  Doppler  frequency  as  that  of  the  Doppler  spectra  shown  in 
Figure  10,  then  the  former  data  can  be  considered  as  if  the  data  were  collected  at  the 
latter's  r/h. 

Mathematically,  if  F(f)  is  the  Fourier  transform  of  f(t) ,  then  F(f  -A/)  corresponds  to  the 
Fourier  transform  of  f(t)eJ27[Aft .  Accordingly,  let  xr(n,m )  and  %  (n,m)  denote  clutter  data 
collected  at  ranges  r  and  r0 ,  respectively,  for  n  =  ().■••„¥  - 1  and  m  =  0, •  ■  ■ , M  - 1 ,  and  let 


A/  =  — 
A 


Jl-(h/r0)2  -Vl -(h/r)2 


(9) 


then  the  operation  of 


Xr '(n,m)  =  xr(n,m)ej2,rASmTr  for  n  =  and  m  =  0,---,M-l  (10) 

compensates  the  Doppler  frequency  of  xM,m)  as  if  x,'(n,m)  were  collected  at  range  r0 .  It 
is  noted  that  the  operation  of  (10)  does  not  alter  spatial  spectra  of  the  data.  Equation  (10) 
may  be  written  as  a  compact  vector  form  of, 

lf=Ar°  q,  (n) 
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where  snapshot  is  an  NMx  1  vector  and  formed  by  stacking  the  data  %r(n,m), 
n  =  0,  •  •  ■ ,  N  - 1 ,  and  M  =().••■,  M  - 1 ;  the  symbol  o  denotes  element-by-element 
multiplication,  and  q(.  the  corresponding  Doppler  compensation  vector  as, 

q,  =sc  ®l  Vxl  (12) 

where  <8>  denotes  Kronecker  product  and, 

sc=[i  exp(/2^A/Tr)  ■■■  exp(/2^(M-l)A/Tr)f  (13) 

In  application,  one  may  designate  r0  =  h ,  the  Doppler  frequency  to  be  compensated 
becomes. 


A/  =  -  —r~  -y/l  -  (A  /  r)  2 
A 


(14) 


The  advantage  of  using  (14)  for  compensating  Doppler  for  range-dependent  forward- 
looking  clutter  is  that  the  Doppler  of  the  mainlobe  clutter  becomes  zero  after  the 
compensation,  which  is  the  same  as  that  of  clutter  collected  from  the  side-looking  mode. 
This  means  all  surface  stationary  targets  in  the  mainlobe  direction  also  have  zero  Doppler 
after  the  compensation.  As  a  demonstration,  the  Doppler  of  the  data  whose  original 
Doppler  spectra  are  shown  in  Figure  10  is  compensated  using  (14)  and  the  resultant 
Doppler  spectra  are  shown  in  Figure  11.  It  can  be  seen  that  after  the  Doppler  compensation 
the  data  appear  as  if  they  were  collected  from  the  side-looking  geometry. 
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Figure  11:  Doppler  frequency  of  the  mainlobe  clutter  shoivn  in  Figure  10  is  shifted  to  zero  after  the 
Doppler  compensation. 
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Therefore,  for  a  non-side-looking  three-dimensional  NxMxK  data  cube,  where  K  is  the 
number  of  range  bins,  we  first  perform  the  Doppler  compensation  from  range  bin  to  range 
bin.  After  the  Doppler  compensation,  the  data  cube  can  be  treated  as  if  it  were  collected 
from  the  side-looking  mode,  then  the  remaining  detection  processing  is  no  different  from 
the  processing  of  the  side-looking  data,  and  any  of  the  fully-adaptive  or  partially-adaptive 
STAP  algorithms  can  be  applied.  The  comparison  between  the  processing  of  side-looking 
data  and  non-side-looking  data  is  shown  in  Figure  12. 


Side-looking 
data  cube 


Non-side-looking 
data  cube 


Output 


Output 


Figure  12:  Non-side-looking  data  cube  will  be  treated  as  side-looking  data  cube  in  STAP  after 
Doppler  compensation. 


The  amount  of  Doppler  compensation  given  in  (14)  is  for  the  forward-looking  case  (0° 
from  nose).  A  more  general  form  of  (14)  is. 


A/  =  -(h/r)2  cos0o 


(15) 


where  </>0  is  the  azimuth  angle  the  array  steers  away  from  the  nose  (-180°  <  (j)()  <180°).  The 
case  of  <f>0  =  ±90°  is  the  side-looking  case  where  no  Doppler  compensation  is  needed. 

Since  all  parameters  of  (15)  are  known  a  priori  for  the  data  to  be  processed,  the 
compensation  using  (15)  should  lead  to  satisfactory  results.  The  effect  of  crabbing  angle  is 
not  considered  in  this  compensation  scheme  and  may  be  treated  /  compensated  separately 
if  necessary.  Kreyenkamp  and  Klemm  (2001)  argued  that  parameters  of  (15)  may  not  be 
known  precisely  and  there  may  be  a  presence  of  crabbing  angle.  They  proposed  a  least 
squares  method  to  estimate  the  Doppler  to  be  compensated.  Unfortunately,  the  proposed 
method  did  not  use  any  radar  data,  so  the  estimate  was  certainly  not  a  least  squares 
estimate.  It  merely  used  a  Taylor  series  expansion  approximation  to  replace  a  precise 
expression. 
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3.  A  Numerical  Example 

A  data  cube  of  sea  clutter  received  by  a  linear  array  operated  in  the  forward-looking  mode 
is  simulated  by  GPARM  (stands  for  Generic  Phased  Array  Radar  Model)  developed  by 
DSTO  (Berry  et  al.  to  be  published).  Assumptions  used  in  the  simulation  include, 

•  Phased  array  is  perfectly  calibrated  and  all  radar  and  platform  parameters  are 
known  (see  Table  1); 

•  NECAPS  model  for  sea  clutter  simulation  (Branson  2000;  Branson  et  al.  2011); 

•  Range-ambiguous  returns  are  included; 

•  Propagation  effects  are  include  (Coleman  and  Berry  2013;  Coleman  et  al.  2009); 

•  Receive  thermal  noise  is  include; 

•  Tapering  windows  to  lower  sidelobes  on  receive  are  used  and  no  tapering  window 
applied  to  transmit. 

The  size  of  the  generated  data  cube  is  NxMxK ,  where  K  =  2433  with  a  range  sampling 
interval  of  15  m  (the  range  resolution  is  30  m),  the  first  range  bin  corresponds  to  a  range- 
to-height  ratio  of  2.0  and  the  last  bin  to  a  ratio  of  9.3. 

The  pulse-based  spatial  spectra  as  well  as  the  channel-based  Doppler  spectra  shown  in 
Figure  13  and  Figure  14  are  for  a  near  range  bin  corresponding  to  a  range-to-height  ratio 
of  2.0.  In  contrast,  the  spatial  spectra  as  well  as  the  Doppler  spectra  shown  in  Figure  15 
and  Figure  16  are  for  a  far  range  bin  corresponding  to  a  ratio  of  9.3.  It  can  be  seen  that 
these  spectra  are  very  similar  to  those  shown  in  Figure  8  to  Figure  11,  although  the 
simulation  used  in  this  section  is  more  sophisticated  and  many  mechanisms  such  as  range- 
ambiguous  returns,  propagation  effects  and  a  more  realistic  clutter  model  are  included. 


Figure  13:  Pulse-based  M  spatial  clutter  spectra  beamfomied  from  the  simulated  clutter  data  for  a 
radar  operating  in  the  forward-looking  mode  with  r/h  =  2.0  (simulated  by  GPARM). 
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Figure  14:  Channel-based  N  Doppler  spectra  of  the  simulated  clutter  data  for  a  radar  operating  in 
the  forward-looking  mode  with  r!h  =  2.0  (Simulated  by  GP ARM) . 
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Figure  15:  Pulse-based  M  spatial  clutter  spectra  beamformed  from  the  simulated  clutter  data  for  a 
radar  operating  in  the  forward-looking  mode  with  r!h  =  9.3  (Simulated  by  GPARM). 


14 


UNCLASSIFIED 


UNCLASSIFIED 


DST-Group-TN-1461 


-30 


-1600  -1200  -800  -400  0  400  800  1200  1600 

Doppler  frequency  (Hz) 

Figure  16:  Channel-based  N  Doppler  spectra  of  the  simulated  clutter  data  for  a  radar  operating  the 
in  forward-looking  mode  with  r/h  =  93  (Simulated  by  GP ARM). 


To  investigate  the  improvement  of  the  Doppler  compensation  to  the  detection 
performance,  the  signal-to-interference-and-noise  ratio  improvement  factor  (SINR  IF)  is 
calculated  for  the  case  of  forward-looking  mode,  using  the  data  without  the  Doppler 
compensation  as  well  as  the  data  after  the  Doppler  compensation. 

The  covariance  matrix  of  clutter  is  estimated  using  the  data  cube  by 

R  =  ^ilklk  (16) 

A  k=l 

The  optimal  weighting  vector  is  (Ward  1994), 

«V=R~‘s  (17) 


where  s  =  s,  ®ss  is  the  space-time  steering  vector;  s,  and  s  v  are  the  temporal  and  spatial 
steering  vectors,  respectively,  and  given  by, 

s,=[l  exp (j2rf,Tr)  ■■■  exp(j2xf ,Tr (M  - l)]r  -fPRF/2<ft  <fPRF/2  (18) 

Ss  =  1  ,vxi  f°r  looking  in  the  mainlobe  direction  (19) 

The  SINR  of  the  optimal  processor  (17)  is  given  by  (Ward  1994), 
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(20) 


where  is  the  single-element  single-pulse  SNR  of  the  target  (i.e.,  the  input  target  signal- 
to-noise  ratio,  SNRln  -  a2  2,  /a1  =  £,). 

The  SINR  IF  of  a  processor  is  defined  as  the  ratio  of  output  SINR  to  input  SNR.  For  the 
above  optimal  processor,  the  SINR  improvement  factor  is. 


SINRif 


o-2|w",s|~ 

w"  Rw" 

opt opt 


(21) 


The  SINR  IF  over  the  whole  Doppler  space  is  calculated  using  the  above  simulated  clutter 
data  cube.  For  comparison.  Figure  17  shows  two  curves  of  SINR  IF,  one  obtained  by 
directly  using  the  data  cube  to  form  the  covariance  matrix  (i.e.  without  the  Doppler 
compensation),  the  other  involving  the  Doppler  compensation.  The  locations  of  the  two 
mainlobe  clutter  notches  are  different,  because  the  one  with  the  Doppler  compensation  has 
made  a  deliberate  Doppler  shift  so  that  the  mainlobe  clutter  has  a  zero  Doppler.  It  can  be 
seen  that  both  curves  reach  the  maximum  SINR  IF  of  10loglo(ATf)  for  frequencies  differing 
from  that  of  the  mainlobe  clutter.  The  important  difference  is  the  width  and  the  depth  of 
the  notch.  For  the  dataset  that  undergoes  the  Doppler  compensation,  the  notch  is  narrower 
and  deeper  providing  better  peak  cancellation  whilst  limiting  cancellation  in  the  adjacent 
frequencies.  This  is  because  that  the  mainlobe  clutter  all  has  a  zero  Doppler  irrespective  of 
range  bins,  so  data  samples  can  be  considered  as  i.i.d.,  and  the  resulting  optimal  weighting 
vectors  can  sufficiently  notch  the  mainlobe  clutter  while  maintaining  the  minimal 
detection  loss  in  the  vicinities  of  mainlobe  Doppler.  On  the  other  hand,  the  data  without 
the  Doppler  compensation  has  varying  Doppler  for  mainlobe  clutter,  and  range  data  are 
not  i.i.d.  As  a  consequence,  the  resulting  weighting  vectors  are  compromised  to  be  the  best 
for  the  overall  dataset,  but  not  the  optimal  from  range  bin  to  range  bin.  Hence  for  near 
range  bins,  detection  loss  occurs  in  the  higher  Doppler  region  whereas  the  detection  loss 
happens  in  the  lower  Doppler  region  for  far  range  bins.  Therefore,  the  notch  for  the  data 
without  the  Doppler  compensation  is  wider.  In  addition,  the  depth  of  the  notch  is  not  as 
deep  as  it  should  be,  and  excessive  false-alarms  will  occur  as  the  processor  cannot 
sufficiently  suppress  the  mainlobe  clutter.  Note  that  even  after  the  Doppler  compensation 
the  notch  is  still  relatively  wide;  this  is  because,  according  to  the  simulation,  it  is  a  case  of 
sea  clutter  that  itself  has  a  Doppler  spectrum  centred  at  zero  Hertz  (i.e.  there  is  internal 
clutter  motion  that  spreads  the  spectrum  of  the  clutter). 
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Figure  17:  S1NR  improvement  factors  with  and  without  Doppler  compensation  (for  the  GPARM 
simulated  dataset). 


The  smaller  the  range  to  height  ratio,  the  more  severe  the  problem,  if  forward-looking  data 
are  directly  processed  using  STAP  algorithms  without  the  Doppler  compensation.  The 
problem,  however,  will  mitigate  for  higher  ratios  of  range  to  height.  According  to  Figure  3, 
once  the  ratio  approaches  10  or  higher,  the  mainlobe  clutter  Doppler  can  be  considered  to 
be  the  same  for  all  range  bins.  In  such  a  scenario  the  Doppler  compensation  in  the  STAP 
may  be  omitted. 

Although  only  a  generic  AEW  scenario  in  a  homogenous  clutter  environment  has  been 
assumed  in  the  above  simulation,  the  described  Doppler  compensation  remains  the  same 
for  heterogeneous  clutter  environments.  Once  the  data  have  been  compensated,  the  data 
can  be  considered  to  be  collected  from  the  side-looking  geometry. 


4.  Range  Ambiguous  Clutter  and  Its  Suppression 

For  airborne  radar  operating  in  medium  to  high  PRF  modes,  range  ambiguity  is 
unavoidable.  Therefore  there  is  a  need  to  address  this  issue.  The  ambiguous  range  is  given 

by. 


ra(k) 


k  =  1,2,— 


(22) 
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where  c  is  the  speed  of  light  and  k  denotes  the  k,h  ambiguous  range,  r  is  the 
unambiguous  range,  and  Tr  =M fPRF  is  the  PRI.  According  to  (4),  if  the  constant 
backscattering  coefficient  cr0  clutter  model  is  assumed,  the  clutter  intensity  is  inversely 
proportional  to  the  range  cubed  (since  the  clutter  patch  area  is  proportional  to  range). 
Alternatively,  if  the  constant  gamma  y()  clutter  model  is  assumed,  because 
cr o  =y0sin 0  =  yoH/r,  where  H  is  the  platform  height,  and  6  is  the  grazing  angle  (the 
grazing  angle  is  approximately  equal  to  the  depression  angle),  the  clutter  intensity  will  be 
inversely  proportional  to  the  range  to  the  4th  power.  As  a  result,  the  propagation 
attenuation  of  the  k"1  ambiguous  range  clutter,  normalised  by  the  unambiguous  range 
clutter,  is. 


A,  =  4  tog 


10 


'“(k) 

V  r  J 


=  4  log 


10 


k  cT 
1  +  A-A4 
2  r 


k  =  1,2,- 


where 


[30  constant  sigma  cr0  clutter  model 
[40  constant  gamma  yt]  clutter  model 


(23) 


Since  the  unambiguous  range  r  is  bounded  by  0  <  r  <  rmax  =  cTr  /  2 ,  the  propagation 
attenuation  Ak  is  a  function  of  unambiguous  range  r  .  Choosing  different  values  of  r ,  we 

can  calculate  values  of  propagation  attenuation  of  the  k,h  ambiguous  range  clutter  which 
are  given  in  Table  2. 

Values  given  in  Table  2,  assuming  clutter  intensity  attenuation  in  a  fashion  of  ( ra(lc)/r )3  or 
(ra(k)  /r)4  may  be  conservative  for  small  grazing  angles.  According  to  experimental 
measurements,  because  of  the  multipath  propagation  effect,  radar  cross-section  (RCS)  of  a 
point  target  generally  attenuates  from  (ra(k)  I  r)A  to  (falk)  l  rf ,  so  that  the  clutter  intensity 

attenuation  can  vary  from  (raik)  / r  f  to  (ra(k)/r)7  for  the  constant  sigma  clutter  model  or 
from  (ra(k)  / r)4  to  (ra(k)  / r)8  for  the  constant  gamma  clutter  model  (Long  2001,  Chapter  4). 


Table  2:  Propagation  attenuation  of  kth ambiguous  range  clutter  (dB),  assuming  clutter  intensity 
inversely  proportional  to  the  range  cubed /  the  range  to  the  4th  poioer). 


Ambiguous  range 

r  =  'Lax  / 1° 

r  —  r  /  2 

max 

r  =  Nax 

k  =  1 

31.2/41.7 

14.3/19.1 

9.0/12.0 

k  =  2 

39.7/52.9 

21.0/28.0 

14.3/19.1 

k  =  3 

44.7/59.7 

25.4/33.8 

18.1/24.1 

k  =  4 

48.4/64.5 

28.6/38.2 

21.0/28.0 
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The  above  discussions  assume  all  range-unambiguous  and  range-ambiguous  clutter  to  fall 
in  the  mainlobe  (or  in  the  same  sidelobe).  In  fact,  situations  can  be  more  complex  if  range- 
unambiguous  clutter  and  range-ambiguous  clutter  fall  in  different  beam  lobes.  For 
instance,  as  depicted  in  Figure  18,  the  range-unambiguous  clutter  falls  in  the  sidelobe 
whereas  the  range-ambiguous  clutter  falls  in  the  mainlobe  or  vice  versa.  It  is  possible  that 
the  range-ambiguous  clutter  can  be  stronger  than  the  range-unambiguous  clutter  when 
radar  is  operated  in  high  PRF  modes. 

Since  Doppler  shifts  of  non-side-looking  clutter  are  range  dependent,  and  the  Doppler 
compensation  scheme  discussed  earlier  only  compensates  the  Doppler  shifts  of  the  range- 
unambiguous  clutter.  Hence,  once  the  range-ambiguous  clutter  is  involved,  the 
suppression  of  range-ambiguous  clutter  becomes  more  crucial  for  radar  operated  in  the 
non-side-looking  mode. 


4.1  Low  Ratio  of  Range  to  Height 

In  general  the  SINR  IF  curve  may  have  two  or  more  notches  when  range-ambiguous 
clutter  is  involved.  Taking  two  notches  as  an  example,  the  second  often  smaller  notch 
corresponds  to  the  weaker  1st  range-ambiguous  clutter,  and  may  locate  at  a  different 
Doppler  depending  on  the  radar's  looking  direction.  For  a  radar  operating  in  the  side¬ 
looking  geometry,  the  Doppler  is  range  invariant,  so  the  second  notch  falls  in  the  first  and 
there  will  be  only  one  notch.  For  a  radar  operating  in  non-side-looking  geometries, 
Doppler  is  range  dependent,  so  the  Doppler  frequencies  of  range-unambiguous  mainlobe 
clutter  and  range-ambiguous  mainlobe  clutter  are  different,  resulting  in  two  or  more 
notches.  If  a  radar  looks  at  near  range  with  a  low  ratio  of  range  to  height,  the  depression 
angles  of  the  ambiguous  range  and  unambiguous  range  are  very  different,  and  one  may 
use  a  common  beamforming  technique  in  the  elevation  to  null  the  illumination  of  the 
range-ambiguous  surface  for  suppression  of  range-ambiguous  clutter.  This  requires  that 
the  transmit  antenna  has  to  be  a  two-dimensional  planar  antenna  so  nulls  can  be  placed  in 
the  required  elevation  beam  pattern. 

For  instance,  for  the  previous  numerical  example,  suppose  radar's  mainlobe  is  steered 
down  at  r  =  1 0  km  with  h  =  5  km  and  r  /  h  =  2 .  According  to  the  given  parameters,  its  first 
unambiguous  range  is  rfl(1)=  56.875 km.  The  difference  of  their  depression  angles  is  25.0° 
which  means  that  the  range-ambiguous  clutter  will  have  to  enter  through  an  elevation 
angle  of  25.0°.  Since  their  elevation  angles  are  very  different,  one  may  utilise  a 
beamforming  technique  to  null  the  range-ambiguous  clutter  to  remove  the  second  notch 
from  the  SINR  IF  curve,  if  the  array  is  a  2D  planar  array.  However,  this  may  be  hard  to  do 
in  practice,  as  there  are  quite  a  few  parameters  determining  the  ambiguous  range  and  the 
elevation  angle  corresponding  to  range-ambiguous  clutter.  Some  radar  systems  may  have 
limited  degrees  of  freedom  in  elevation  for  beamforming. 
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4.2  High  Ratio  of  Range  to  Height 

If  a  radar  steers  at  far  range  with  a  high  ratio  of  range  to  height,  the  difference  of 
depression  angles  of  unambiguous  range  and  ambiguous  range  is  small,  and  it  is  perhaps 
difficult  to  null  the  range-ambiguous  clutter  using  the  beamforming  technique.  However, 
at  a  high  ratio  of  range  to  height,  the  Doppler  of  the  clutter  becomes  approximately 
independent  of  range,  and  the  range-ambiguous  clutter  has  the  same  Doppler  as  that  of 
the  range-unambiguous  clutter,  similar  to  the  side-looking  case.  Hence  the  detection  loss 
caused  by  the  range-ambiguous  clutter  is  negligible.  The  only  difference  is  that  the 
involvement  of  the  ambiguous  clutter  increases  the  rank  of  the  clutter  covariance  matrix, 
so  it  requires  more  sample  data  to  estimate  the  covariance  matrix  in  order  to  achieve  the 
desired  performance. 


4.3  Medium  Ratio  of  Range  to  Height 

It  is  difficult  to  null  the  range-ambiguous  clutter  when  a  radar  is  operated  at  a  medium 
ratio  of  range  to  height,  as  the  difference  in  depression  angles  is  small.  However,  the 
difference  of  the  corresponding  Doppler  frequencies  is  also  small.  Overall  the  involvement 
of  ambiguous  clutter  slightly  broadens  the  clutter  notch,  causing  some  minor  detection 
loss  in  the  vicinities  of  the  Doppler  of  mainlobe  clutter,  which  we  may  have  to  bear  with 
unless  other  techniques,  such  as  suppression  of  range-ambiguous  clutter,  are  further  used. 


5.  Potential  Techniques  for  Suppression  of  Range- 

Ambiguous  Clutter 


5.1  3D  Data  Processing 

Commonly  STAP  algorithms  perform  two-dimensional  space-time  adaptive  processing. 
The  space  involved  normally  refers  to  the  azimuth  direction.  In  fact,  airborne  phased 
arrays  are  typically  of  a  rectangular  latticed  structure.  Data  from  antenna  elements  in 
elevation,  if  available,  provide  another  dimension  of  processing  leading  to  the  so-called  3D 
STAP.  Adding  elevation  processing  under  some  conditions  can  offer  significant  clutter 
suppression  improvement,  allowing  further  suppression  of  the  interference  sources  that 
have  identical  azimuth  and  Doppler  but  different  elevation  from  the  target  signals  (Fertig 
and  Krich  2005;  Dong  2011). 

The  achievement  of  further  clutter  suppression  by  adding  elevation  processing  is  twofold. 
First,  for  close  range  bins,  airborne  targets  and  the  ground  clutter  may  have  significant 
elevation  angles  as  shown  in  Figure  18.  Adaptive  beamforming  processing  in  elevation  can 
thus  separate  target  signals  from  undesired  interferences,  proving  clutter  suppression  in 
the  elevation  domain.  Secondly,  when  radar  is  operated  in  medium  to  high  PRF  modes, 
the  range  is  often  ambiguous  as  mentioned  earlier.  In  some  cases  as  discussed  previously, 
nothing  can  be  done  to  mitigate  the  effect  of  the  range-ambiguous  clutter.  However,  since 
the  range-ambiguous  clutter  comes  from  a  different  elevation  angle,  a  proper  adaptive 
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processing  in  the  elevation  domain  can  filter  out  the  ambiguous  range  clutter  to  improve 
SINR  in  Doppler  regions  that  are  otherwise  degraded  by  the  range-ambiguous  clutter.  In 
their  paper,  Fertig  and  Krich  (2005)  show  an  example  in  which  the  range-ambiguous 
clutter  is  suppressed  and  target  detection  improved  by  3D  STAP. 

In  principle,  2D  STAP  algorithms  can  be  readily  extended  to  3D  without  any  difficulty 
(Fertig  and  Krich  2005).  The  only  problem  is  that  adding  another  dimension  in  STAP 
significantly  increases  the  scale  of  data  recording  and  processing,  and  enlarges  the  size  of 
the  covariance  matrix.  Flence  the  requirements  for  sample  data  support  and  computational 
budget  also  significantly  increase  accordingly. 


Figure  18:  Elevation  angles  of  an  airborne  target  and  the  ground  clutter  can  be  significantly 
different  for  close  range  targets. 


5.2  Frequency  Diversity  Technique 

The  aforementioned  3D  STAP  (adaptive  in  azimuth,  elevation  and  Doppler)  greatly 
increases  the  computational  cost  and  demands  much  more  sample  data  to  support.  In 
practice,  there  may  not  be  sufficient  i.i.d.  samples  for  estimating  the  covariance  matrix2. 
Recently  a  concept  for  using  frequency  diversity  to  mitigate/ suppress  range-angle 
dependent  clutter  has  been  proposed  (Baizert  et  al.  2006;  Antonik  et  al.  2006b;  Sammartino 
et  al.  2013).  Frequency  diverse  array  (FDA)  differs  from  its  conventional  phased  array 
counterpart  by  employing  a  small  frequency  increment  across  array  elements.  It  can 
generate  a  range,  angle  and  time  dependent  beam  pattern  (Antonik  et  al.  2006b;  Huang  et 
al.  2008;  Wang  2013).  The  small  frequency  increment  provides  additional  degrees  of 
freedom  (DoF),  which  lead  to  a  separation  of  range-ambiguous  clutter  from  range- 
unambiguous  clutter. 


2  It  is  possible  to  use  previous  CPI  data  for  estimating  the  covariance  matrix  for  the  current  CPI  to 
overcome  the  sample  data  issue  (see  Dong,  Y.  (2005),  "Approximate  invariance  of  the  inverse  of  the 
covariance  matrix  and  the  resultant  pre-built  STAP  processor",  DSTO-RR-0291). 


UNCLASSIFIED 


21 


UNCLASSIFIED 

DST-Group-TN-1461 

Using  FDA  to  suppress  the  range-ambiguous  clutter  is  explained  in  the  following.  For  a 
FDA  the  waveform  radiated  from  an  array  element  differs  from  the  others  by  a  carrier 
frequency  increment  A f  (Antonik  et  al.  2006a).  For  a  linear  array, 

./„  =  ./o+«A/'  n  =  0,-,N-l  (24) 


Assume  all  N  waveforms  are  mutually  orthogonal  (de  Wit  et  al.  2011),  and  hence  the 
received  signals  can  be  separated  in  frequency.  In  order  to  re-achieve  the  transmit  gain,  the 
received  signal  is  synthesised  prior  to  further  processing.  In  other  words,  transmit  is 
realised  on  an  element  basis  whilst  receive  is  achieved  using  the  whole  aperture.  Signals  of 
different  frequencies  are  then  separated. 

Let  us  discuss  the  spatial  frequency  of  FDA.  Consider  that  the  geometry  between  the  array 
and  a  distant  scatterer  is  defined  by  range  r  (referenced  to  channel  0),  azimuth  angle  </> 
and  depression  angle  0.  The  received  signal  by  the  full  aperture  for  the  waveform 
transmitted  by  the  nth  channel  is  (common  terms  are  omitted). 


Yn  (&’  fc  r)  =  eXP 


,2  nfn 


(r  +  Wcos<9sin  </>) 


N- 1 

Z  exp 

k= 0 


j  — ^-(r  +  k  d  cos  6  sin  (p) 


n  =  0,---,N -1  (25) 


where  the  first  term  is  the  phase  corresponding  the  waveform  emitted  from  the  nth 
channel  to  the  scatterer,  and  the  second  sum  term  is  the  returned  waveform  synthesised  by 
the  whole  aperture.  The  summation  gives. 


Yn(0,<l>-,r)  =  exp 


,2  nfn 


2 /-  +  ]  — — -  +  /7  \d  cos  9  sin  (f> 


-i  sin 


Nn  —  dcosOs'mtj) 
c 


I  f 

J  sin  n— Ld  cosOsmcp 

V  c 


! =  0,  *  ■  ■ ,  N  —  1  (26) 


Suppose  the  radar  steers  to  <p  =  0 ,  the  beamforming  of  N  channels  gives. 


n- i  n- i  (  2  nf, 

m0 ;/■)=  Zn.(^0;r)  =  iVZ  exp  j  ~^2r 


k= 0 


k= 0 


=  A^exp 


,2^/o 


.2  jd±f 


j - -2r  exp  j - —{N  —  l)r 


•  i  2M/ 
sin  N - r 


.  (2  xAf 


sin 


V  c 


(27) 


Obviously,  for  A/  =  0 ,  (27)  simplifies  to  N2  exp 


.2^o 


2  r 


,  an  expected  result  for  the 


conventional  phased  array,  where  N2  is  the  round  trip  coherent  processing  gain  (i.e.,  a 
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transmit  gain  of  N  and  receive  gain 
delay. 


of  N),  and  exp 


As  shown  in  (27),  the  beam  pattern  in  range3  |r(0,O;r)|  periodically  reaches  maximum 
when  sin(2^r  A/ / c)  =  0  .  Locations  of  |  V(9.()\r)  |  peaks  are. 


r  = 


knc 
2  nAf 


kc 

2Af 


k  =  0,  I,---, 


(28) 


Accordingly,  |r(0,O;r)|  periodically  reaches  maxima  at  intervals  of  Ar  =  c/(2Af).  Figure  19 
shows  the  beam  pattern  in  range  periodically  varying  with  range.  Note  that  the  range 
shown  in  the  figure  is  range  relative  to  the  range  of  interest  (beam  patterns  in  range  are 
discussed  in  the  far  field  region).  For  instance,  if  we  are  interested  in  the  radar 


measurement  at  r  =  r0 ,  a  phase  term  of  exp 


Mn 


■J- 


2rn 


should  be  compensated  for 


v  c 

yn(9,<p-,r) ,  n  =  0,- •, N - 1 ,  so  that  the  radar  wave  will  be  focused  at  r  =  r0 ,  and  the  associated 
beam  pattern  in  range  will  be  the  same  as  what  is  shown  in  Figure  19,  and  the  relative 
range  0  refers  to  r  =  r0 . 


Figure  19:  Beam  -pattern  in  range  (N  =  29). 


3  The  beam  pattern  of  traditional  radar  normally  refers  to  the  beam  pattern  in  azimuth  as  its  beam 
pattern  is  independent  of  range.  The  beam  pattern  of  FDA  is  however  dependent  on  both  azimuth 
and  range.  The  beam  pattern  |  Y(9S)\r)  \  discussed  here  is  the  beam  pattern  in  range,  and  hence  the 
term  "beam  pattern  in  range"  is  used. 
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It  can  be  seen  that  compared  to  a  conventional  phased  array  radar  whose  beam  pattern  in 
range  is  independent  of  range,  the  beam  pattern  in  range  of  a  FDA  radar  varies 
periodically  in  range  adding  another  dimension  of  DoF.  One  of  its  potential  applications  is 
the  suppression  of  the  range-ambiguous  clutter.  We  recall  that  the  ambiguous  range 
interval  is  cTr!2  (see  (22)).  In  order  to  focus  at  the  unambiguous  range  r ,  a  phaser  of 


exp 


;2  Tfn 


Il¬ 


ls  needed  to  compensate  yn,  «  =  0, •  •  ■ , N - 1 .  After  the  compensation,  the 


relative  range  between  the  focused  range  and  the  first  ambiguous  range  is  cTr/2 .  Because 
our  goal  is  to  suppress  the  range-ambiguous  clutter,  we  may  choose  such  parameters,  so 
that  the  first  ambiguous  range  falls  in  the  middle  of  the  two  maxima  of  the  beam  pattern  of 
Figure  19.  This  criterion  requires. 


c 

2A J 


(k,.+ 1/2) 


(29) 


where  kr  is  a  natural  number  of  kr  =  0,1,  ■  ■  • .  From  (29)  one  obtains. 


fr 


A/ 

kr  +  0.5 


(30) 


Accordingly  one  selects  an  appropriate  integer  kr  in  such  a  way  so  that  the  resulted  fPRF 
is  close  to  the  desired  one.  With  such  a  selection,  the  clutter  return  from  the  1st,  3rd,  5th  and 
so  on  ambiguous  range  will  be  suppressed.  Alternatively,  one  may  select  fPRF  by. 


// 


PRF 


A/ 

kr  +  0.25 


(31) 


By  this  way,  the  clutter  returns  from  the  1st,  2nd,  3rd  and  5th,  6th,  7th  ambiguous  range  and  so 
on  will  be  suppressed.  As  an  example,  differences  between  the  selected  fPRF  and  the 
desired  fPRF  based  on  (30)  and  (31)  are  given  in  Table  3. 


Table  3:  Choose  proper  PRFs  to  suppress  range-ambiguous  clutter  (A/  =  100  kHz ). 


Desired  PRF 

Selected  PRF  based  on  (30) 

Selected  PRF  based  on  (31) 

1000  Hz 

995  Hz 

998  Hz 

2000  Hz 

1980  Hz 

1990  Hz 

3000  Hz 

2985  Hz 

3008  Hz 

In  addition  to  the  suppression  of  range-ambiguous  clutter,  there  is  another  potential 
application  of  FDA.  When  radar  is  operated  in  medium  to  high  PRFs,  range  becomes 
ambiguous.  In  the  conventional  phased  array  radar,  whether  a  target  presents  in 


24 


UNCLASSIFIED 


UNCLASSIFIED 

DST-Group-TN-1461 

unambiguous  range  or  ambiguous  range  is  not  distinguishable,  and  the  radar  needs  to  use 
different  PRFs  to  discriminate.  With  the  FDA,  it  requires  a  unique  phase  compensation  to 
focus  at  the  range  of  interest.  As  a  result,  if  the  phase  compensation  is  for  the 
unambiguous  range,  a  target  in  the  ambiguous  range  will  be  greatly  attenuated  or 
suppressed.  To  detect  a  target  located  in  the  ambiguous  range,  the  associated  phase 
compensation  has  to  be  used.  Therefore,  there  is  perhaps  no  need  for  DFA  to  collect 
different  PRF  data  to  determine  the  true  location  of  target  even  though  the  radar  is 
operated  in  a  range-ambiguous  PRF. 

In  the  above,  the  phase  compensation  is  carried  out  in  the  radio  frequency  domain,  which 
may  not  be  convenient.  In  fact,  the  phase  compensation  can  be  carried  out  in  the  baseband. 
The  phase  difference  between  nth  channel  and  the  reference  channel  (0th  channel)  is, 


A  y/n  (0,  </),  r )  =  Z  yn  {0,4,r)-Z  y0  {0,  </>,  r) 


2k  n  A /  _  2k 

- —  2  r  +  — 

C  An 


1  + 


Affjv-l 

/ol  2 


+  n 


n  d  cos  9  sin  <j> 


(32) 


If  there  is  no  frequency  diversity.  A/  =  0 ,  the  first  term  vanishes,  and  the  phase  difference 


2k 


becomes  the  familiar  one  of  — nd cos# sin ^  for  a  conventional  phased-array  radar,  which  is 


independent  of  range.  Equation  (32)  may  be  written  as. 


=  2k  fr  +  2k  fs 


(33) 


Equation  (33)  shows  that  with  a  FDA,  the  spatial  frequency  consists  of  two  components: 
the  first  one  is  range  dependent,  and  the  second  one  is  orientation  dependent.  In  the 

processing,  a  proper  phaser  of  exp|^-  j  2r  j ,  n  =  0,1,  ■  ■  •,  N  - 1 ,  has  to  be  compensated  for 

the  range  of  interest,  so  that  the  waveform  focuses  at  the  range  of  interest,  and  range- 
ambiguous  clutter  is  suppressed.  The  remaining  processing  will  be  no  different  to  the 
conventional  phased-array  data,  rather  than  using  a  different  spatial  steering  vector.  All 
STAP  algorithms  developed  for  the  conventional  phased-array  can  be  applied. 

There  is  a  need  to  justify  that  the  Doppler  properties  of  clutter  is  approximately 
unaffected.  For  a  clutter  patch  located  at  {9,4),  the  associated  Doppler  for  the  72th  channel 
is. 


fan  =  -^-cos^sin^  =  — ^-/o(l  +  wA/7 /0)cos#sin^  «  ^^-cos0sin^ 
A„  c  A0 


(34) 


Therefore,  the  effect  of  the  frequency  diversity  on  the  Doppler  is  negligible. 
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In  summary,  a  block  diagram  of  the  FDA  data  processing  technique  to  suppress  range- 
ambiguous  clutter  is  depicted  in  Figure  20.  The  main  differences  between  the  FDA  data 
processing  and  the  conventional  phased-array  data  processing  include. 


•  Array  synthesis  is  performed  on  receive  for  FDA,  whereas  it  is  automatically 
performed  at  transmit  for  a  phased-array; 

•  Phase  compensation  is  required  for  FDA  data  to  focus  the  radar  wave  at  the  range 
of  interest  whereas  no  such  a  step  is  required  for  the  phased-array  data; 

•  The  spatial  steering  vector  is  different. 

The  remaining  STAP  processing  is  identical  for  both  FDA  and  phased-array  data. 


FDA 

data 


Output 


Figure  20:  Block  diagram  of  FDA  data  processing  to  suppress  range-ambiguous  clutter. 


Detailed  architectures,  characteristics  and  techniques  of  the  FDA  radar  will  be  investigated 
in  a  separate  report. 


6.  Conclusions 


This  Technical  Note  is  a  supplement  to  an  earlier  DSTO  Technical  Note,  "Overview  of 
STAP  algorithms"  (Dong  2011).  It  is  aimed  at  the  processing  of  phased  array  data  collected 
from  non-side-looking  geometries.  In  general,  STAP  algorithms  developed  for  a  radar 
operating  in  the  side-looking  mode  cannot  be  directly  applied  to  data  collected  by  a  radar 
operating  in  non-side-looking  modes,  including  the  forward-looking  mode.  A  simple 
Doppler  compensation  scheme  is  proposed,  so  that  the  Doppler  of  mainlobe  clutter  is 
shifted  to  zero  and  becomes  range  invariant,  like  the  clutter  data  collected  from  the  side¬ 
looking  geometry.  After  the  Doppler  compensation,  all  STAP  algorithms  developed  for  the 
side-looking  geometry  can  be  applied. 

Also  briefly  discussed  are  a  few  potential  techniques  for  the  suppression  of  range- 
ambiguous  clutter,  including  beamforming  techniques  and  the  3D  STAP.  The  frequency 
diversity  technique  discussed  in  literature  is  interesting  but  further  study  is  needed.  A 
detailed  study  of  frequency  diverse  array  will  be  reported  separately  in  the  future. 
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